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Numerical studies of confinement in the Landau gauge 
by the larger lattice simulation 
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The infrared behaviour of QCD in the Landau gauge is studied in 24 4 and 32 4 lattices by measuring the 
Kugo-Ojima parameter, gluon propagator, ghost propagator and the QCD running coupling. The Kugo-Ojima 
parameter c remains about 0.8 instead of the expectation value 1. We observe infrared suppression of the singu- 
larity of the gluon propagator and enhancement of the ghost propagator but the shifts are smaller than those of 
Dyson-Schwinger approach. The momentum dependence of the running coupling in the infrared region is weak 
and Q s (0) ~ 1. 



1. The infrared behaviour of QCD in the 
Landau gauge 

The colour confinement in QCD can be formu- 
lated as an absence of coloured asymptotic states 
in the BRST formalism, and the presence of long- 
range 1/p 4 correlation between colour sources. In 
the Gribov-Zwanziger's theory of Landau gauge 
QCD[1], the long-range correlation is attributed 
to the ghost propagator, i.e. the inverse of the 
Faddeev-Popov(FP) matrix. The ghost propaga- 
tor is more singular and the gluon propagator is 
less singular than the tree level in the infrared 
region. 

In Kugo-Ojima's theory [2], a sufficient condi- 
tion for an absence of coloured asymptotic states 
is that the parameter u ab defined by the two-point 
function of the FP ghost fields, c(x),c(y), and 
A v (y)[3], satisfies u ab (0) = -S ab . 

2. The Kugo-Ojima confinement criterion 
and the Gribov-Zwanziger's theory 



Let two-point tensor G^ xy ab be 
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ab 
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and the horizon function H(U) 

H(U) = J2 G^v aa - (N 2 - WW) (2) 

x,y,a 

where N is the number of colours, now N = 3. 
E(U) depends on the options of the definition of 
gauge field, [/-linear or log U type; 
[/-linear: E(U) = J2i w Re trU ^ 
logU: E(U) = £ ; , a tx (\^S(Ai)\ a ) /(N 2 -l), 
where Ai = adjA t , and S(x) 
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th(x/2) ' 

One can show that in the core region [1], the 
horizon function H{U) is negative and in the in- 
finite volume limit limy^oo (H(U))/V = which 
is called horizon condition. 

Zwanziger showed in his lattice formulation 
that the ghost self-energy cancels precisely the 
zero-th order p 2 term and the ghost propagator 
behaves as 1 /p 4 in the infrared asymptotic region 
provided no additional non-perturbative effect is 
involved [1]. Recent results of Dyson-Schwinger 
approach suggest the nonperturbative aspects are 
in effect. 

Kugo also showed in the continuum formula- 
tion that the ghost propagator G(p 2 )/p 2 behaves 
in the infrared asymptotic region as [2] 

G(p 2 )/p 2 ~l/{p 2 (l + u(0) + O(p 2 ))} (3) 

where u(p 2 ) ab — S ab u(p 2 ). Thus Kugo-Ojima 
condition leads to the same singular behavior as 
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Table 1 

The Kugo-Ojima parameter c , trace e/d and the horizon function deviation factor h in U— linear(suffix 
1) and logt/(sufhx 2) version, (3 = 6.0. The bottom line is of (3 = 6.4. 



L 


Cl 


e x /d 


K 


C2 


e 2 /d 


h 2 


16 


0.576(79) 


0.860(1) 


-0.28 


0.628(94) 


0.943(1) 


-0.32 


24 


0.695(63) 


0.861(1) 


-0.17 


0.774(76) 


0.944(1) 


-0.17 


32 


0.706(39) 


0.862(1) 


-0.15 


0.777(46) 


0.944(1) 


-0.16 


32 


0.650(39) 


0.883(1) 


-0.23 


0.700(42) 


0.953(1) 


-0.25 



Zwanziger's, and derives the universal renormal- 
Z 1 

izaiton factors as — = — = 1 + u(Q) = 0. 
Z 3 Z 3 

Putting Kugo-Ojima parameter as u(0) = —c, 
one finds that the horizon condition is written as 
(£)+(d-l)c-e = (d-1) (c-j) = (d-l)h = 0. 

where e = E/V and d = 4. 

Thus one sees that both the above arguments 
appear to be consistent with each other provided 
the lattice covariant derivative meets with the 
continuum one in the naive sense, i.e., e/d = 1. 

Kugo discussed in his lecture [2] that the colour 
confinement can be realized via either (1) Z\ = 
and Z 3 =finite or (2) Z x /Z 3 = but/and Z 3 = 

Zwanziger[l,4] argues that the second possibil- 
ity is realized in the SU(2) colour case. 

Our numerical results are summarized in Ta- 
ble 1. We obtained that u£(0) is consistent to 
— cSi, c = 0.8 in SU(3) quenched simulation, 
(3 = 6.0, on 24 4 and 32 4 . 

In the case of f3 = 6.4 on 32 , the value of 
e/d becomes closer to 1, but \h\ becomes larger 
as compared to those of (3 = 6.0, and thus c be- 
comes smaller. We need to perform larger lattice 
simulation to approach the continuum limit. 

3. The two-point functions (gluon propa- 
gator, ghost propagator) 

We measured gluon propagator of (3 = 6.0 on 
lattice 16 4 , 24 4 and 32 4 , and observed clear ten- 
dency of the infrared suppression in the case of 
24 4 and 32 4 . 

The ghost propagator is infrared divergent and 
its singularity can be parametrised as p~ 2 ( 1+aa > , 

where p 2 = V(4sin 2 ^), (-L/2 < k„ < L/2). 
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Figure 1. The log of the gluon dressing function 
as a function of the log of the lattice momentum 
pa. (3 = 6.0, 32 4 in [/-linear version. 



In the 32 4 lattice, the overall scale difference 
in the log U and U— linear versions disappear in 
gluon propagator, but it still remains in ghost 
propagator. 

The gluon propagator and the ghost propaga- 
tor are evaluated also from Dyson-Schwinger ap- 
proach [5]. They found qualitatively similar solu- 
tions for dressing functions of the gluon and the 
ghost, Z{p 2 ) - {p 2 )- a ° and G{p 2 ) - (p 2 )~ a ^ 
with a.r> + 2ac = 0, qg ~ 0.595 and an — 
-1.19[5]. 

The gluon propagator vanishes at momen- 
tum if ac > 0.5. Our lattice simulation of 
the ghost dressing function G(p 2 ) shows that the 
value Qg ~ 0.18 remains unchanged at 1/3 of the 
Dyson-Schwinger one through 16 4 to 32 4 simu- 
lations. The negative ao appears in 24 4 lattice 
and ctD — —0.4 in 32 4 lattice. We observe the 
tendency that |od| increases as the lattice size 



3 




Figure 2. The log of the ghost dressing function as a function of the log of the lattice momentum pa(left) 
and the running coupling a s as a function of the lattice momentum p(GeV) (right). Triangles are in log U 
and diamonds arc in [/-linear version of (3 = 6.0, 32 4 lattice. 



increases, but our data still stays about 1/3 of 
the prediction of the Dyson-Schwinger approach. 

In contrast to the first data of the ghost propa- 
gator [6], suppression at the lowest momentum is 
not seen in our data of the ghost propagator. For 
calculation of the ghost propagator, we adopted 
the multigrid Poisson solver whose accuracy was 
kept within 10~ 5 , and we set 1% as an ending 
condition of perturbative calculation of the in- 
verse FP operator. The Landau gauge accuracy 
of divA = is 10 in maximum norm. 

Using the ghost and gluon dressing func- 
tion, the running coupling can be evaluated 
from a s (pa) = (g 2 (a)/4:iT)Z(pa)G(pa) 2 ~ 
(pa)~ 2 " D ~ 4aG , where a is the lattice spacing. In 
the case of j3 = 6, 32 4 , which has the smallest \h\, 
the suppression of the gluon dressing function due 
to the loss of correlation of the gluon field and the 
enhancement of the ghost dressing function can- 
cel. We find an + ^olq — near the lowest mo- 
mentum point and a s (0) ~ 1. The infrared sup- 
pression of the running coupling measured from 
the three point function of gluon [7] and from the 
ghost-antighost gluon [8] in MOM scheme would 
be attributed to the loss of correlation of gluon 
fields. 



4. Discussion and outlook 

In the (3 = 6, 32 4 lattice simulation \h\ is about 

0. 15 and we need extrapolation to \h\ = 0. There 
remain finite size effects in ghost propagator in 
32 4 lattice, and also in the Kugo-Ojima param- 
eter which manifest themselves in their depen- 
dence on the definition of the gauge field. The 
continuum limit should not depend on the defini- 
tion of the gauge field and so a systematic method 
for extrapolating the continuum limit of c in the 
two definitions is necessary. 

This work was supported by the KEK super- 
computing project No. 02-82. 
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